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ABSTRACT 

Supernova remnants may exhibit both thermal and nonthermal X-ray emis- 
sion. In a previous study with ASCA data, we found that the middle-aged 
supernova remnant RCW 86 showed evidence for both processes, and predicted 
that observations with much higher spatial resolution would distinguish harder 
X-rays, which we proposed were primarily synchrotron emission, from softer, 
thermal X-rays. Here we describe Chandra observations which amply confirm 
our predictions. Striking differences in the morphology of X-rays below 1 keV 
and above 2 keV point to a different physical origin. Hard X-ray emission is corre- 
lated fairly well with the edges of regions of radio emission, suggesting that these 
are the locations of shock waves at which both short-lived X-ray emitting elec- 
trons, and longer-lived radio-emitting electrons, are accelerated. Soft X-rays are 
spatially well-correlated with optical emission from nonradiative shocks, which 
are almost certainly portions of the outer blast wave. These soft X-rays are well 
fit with simple thermal plane-shock models. Harder X-rays show Fe Ka emis- 
sion and are well described with a similar soft thermal component, but a much 
stronger synchrotron continuum dominating above 2 keV, and a strong Fe Ka 
line. Quantitative analysis of this line and the surrounding continuum shows that 
it cannot be produced by thermal emission from a cosmic-abundance plasma; the 
ionization time is too short, as shown both by the low centroid energy (6.4 keV) 
and the absence of oxygen lines below 1 keV. Instead, a model of a plane shock 
into Fe-rich ejecta, with a synchrotron continuum, provides a natural explana- 
tion. This requires that reverse shocks into ejecta be accelerating electrons to 
energies of order 50 TeV. We show that maximum energies of this order can be 
produced by radiation-limited diffusive shock acceleration at the reverse shocks. 
In an Appendix, we demonstrate that an explanation of the continuum as due to 
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nonthermal bremsstrahlung is unlikely. 

Subject headings: ISM: individual (RCW 86) - supernova remnants - X-rays: 
ISM 

1. NONTHERMAL X-RAYS AND RCW 86 

While most shell supernova remnants (SNRs) show thermal X-ray spectra dominated 
by emission lines of highly ionized species of C, N, O, Ne, Mg, Si, S, and Fe, a few 
show featureless spectra which can be well understood as synchrotron emission: SN 1006 
(Koyama et al. 1995; Dyer et al. 2001); G347.3-05 (Koyama et al. 1997; Slane et al. 1999), 
G266.2-1.2 (Slane et al. 2001), and G28.6-0.1 (Bamba et al. 2001). The absence of lines 
means that a synchrotron interpretation is almost unavoidable; impossibly low abundances 
and/or peculiar physical conditions would be required to suppress lines completely in 
a thermal plasma (e. g., Hamilton, Sarazin, & Szymkowiak 1986; Laming 1998), while 
nonthermal bremsstrahlung implies the presence of electrons that should excite lines as 
well, and inverse-Compton scattering of any photon population would result in far too 
hard a spectrum. If synchrotron emission can dominate the spatially-integrated X-ray 
spectrum of a few remnants, it is likely to play some role in many more, e. g., by totally 
dominating thermal emission in selected locations within them, or by overwhelming thermal 
continua (but not lines) in thermally-emitting regions. The signature of an X-ray spectrum 
with both thermal and synchrotron radiation would be unusually weak lines, diluted by 
the synchrotron continuum. We interpreted weak X-ray lines in the remnant RCW 86 
by such combination of thermal and synchrotron emission, based on ASCA observations 
(Borkowski et al. 2001, hereafter BRRD). A very similar interpretation of RCW 86 spectra 
was also proposed by Bamba, Koyama, & Tomida (2000), based on the ASCA Performance 
Verification data. 
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RCW 86 (G315.4-2.3) is a large (42' in diameter) shell-like SNR with a moderate 
blast wave speed in the range of 400-900 km s -1 , determined from observations of 
Balmer-dominated shocks which almost completely encircle the remnant (Long & Blair 
1990; Smith 1997; Ghavamian et al. 2001). At a kinematic distance of 3 kpc, its large size 
and the moderate blast wave speed imply a mature (~ 10 4 yr old) remnant (Rosado et al. 
1996). The bright complex of optical emission filaments at its SW corner, to which the 
optical designation "RCW86" actually refers, has been generally interpreted as a region 
where the blast wave impacted a dense (~ 10 cm" 3 ) interstellar cloud. This is also where 
X-ray and radio emission are the brightest. Radio observations by the MOST telescope 
at 843 MHz (Whiteoak & Green 1996), with an angular resolution of about 45", and by 
the Australia Telescope Compact Array (ATCA) at 1.3 GHz (Dickel et al. 2001), with 8" 
resolution, show a complete shell of varying brightness. The radio morphology in the SW 
region of the remnant is somewhat unusual, consisting of a bright ridge of emission mostly 
interior to optical and soft X-ray filaments which mark the location of the blast wave. 
ASCA Performance Verification X-ray spectra showed strikingly weak X-ray lines, much 
weaker than expected from a normal (solar) abundance plasma, with the exception of the 
Fe Ka line at 6.4 keV (Vink et al. 1997). These results were confirmed by deeper ASCA 
observations (BRRD) and by more recent BeppoSax observations (Bocchino et al. 2000). 
A substantial Fe Ka equivalent width in the SW suggests enhanced Fe abundance with 
respect to solar, irrespective of the origin of the underlying high-energy continuum, which 
led Bocchino et al. (2000) to consider a possibility that we are observing metal-enriched 
supernova (SN) ejecta. 

In our earlier work (BRRD), we used long ASCA observations to deduce the presence of 
three spectral components: a soft thermal component (kT ~ 0.8 keV), with a relatively long 
ionization timescale (r = n e t ~ 2 x 10 11 s cm" 3 ) and with approximately solar abundances, 
modeled as a plane shock (XSPEC model pshock); a synchrotron component, modeled as 
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the high-energy tail of the radio spectrum with XSPEC model srcut; and a hot thermal 
component (kT ~ 5 keV, r ~ 5 x 10 s s cm -3 ) required to account for Fe Ka emission. 
These components varied in relative strength in three regions in the southwest corner 
of RCW 86. The srcut model is just the synchrotron emissivity of a power- law electron 
distribution with an exponential cutoff, calculated numerically (Reynolds & Keohane 1999); 
it is characterized by a single free parameter, the peak frequency f ro iiofr radiated by electrons 
with the e-folding energy of the exponential cutoff. The fitted values of f ro iiofr varied from 
1 to 4 x 10 16 Hz, implying (for B ~ 10 /iG) maximum electron energies of order 20 TeV. 
The synchrotron fitting used as inputs the radio fluxes for each region measured from the 
MOST image (Whiteoak & Green 1996), and an assumed constant spectral index of 0.6. 

This analysis explained the ASCA data reasonably well, but with ASCA's spatial 
resolution of several arcminutes, we were unable to perform several more rigorous tests of the 
hypothesis of synchrotron X-ray emission. Our spectral fits suggested that the soft thermal 
emission is due to nonradiative shocks as identified by Ha emission, while the synchrotron 
X-rays should show a spatial distribution similar to radio emission. In an attempt to verify 
these predictions, we observed RCW 86 with the Chandra X-ray Observatory in February 
2001. We report here on these observations, which not only confirm our predictions, but 
allow us to associate synchrotron X-ray emission with Fe-rich SN ejecta seen through Fe 
Ka line emission, and identify all high-energy continuum as synchrotron emission. 



2. CHANDRA OBSERVATIONS 

RCW 86 was observed for 92 ks February 1-2, 2001 (observation id 500170) by the 
Chandra X-ray Observatory (Weisskopf, O'Dell, & van Speybroeck 1996) Advanced CCD 
Imaging Spectrometer (ACIS: Garmire et al. 1992). Chips 12, 13, SO, SI, S2 and S3 were 
used, with the pointing center on S3. The focal-plane temperature was —120° C. The 
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observations were corrected for radiation damage (referred to as charge transfer inefficiency, 
CTI) that took place early in Chandra 's mission, following the procedure in Townsley et 
al. (2000) which corrects simultaneously for both position-dependent gain shifts and event 
grade changes. The data were extracted and processed using the Chandra Interactive 
Analysis of Observations software package (CIAO) version 2.1.1. The images were corrected 
with exposure and instrumental efficiency maps to address chip gaps and mirror vignetting. 
Instrumental efficiency maps for each chip were created for each energy band (0.5 - 1.0, 1.0 
- 2.0, and 2.0 - 8.0 keV), for a single energy in the midrange of each band. The exposure 
map has been correctly weighted for each chip. The Chandra images shown here are binned 
by a factor of 4 to a resolution of 2"/pixel. 

We analyzed extracted spectra with XSPEC 11.1.0 (for introduction to XSPEC, 
see Arnaud 1996 1 ). The spectra were grouped by a minimum of 15 counts per bin. We 
have experimented with the source spectra and background, comparing CTI-corrected 
spectra, using Pennsylvania State University software tools (Townsley et al. 2000 2 ), and 
non- CTI-corrected spectra, using standard CXC CIAO tools. The CTI-corrected spectra 
give more accurate centroids of lines such as Mg and Si, resulting in a particularly significant 
improvement for the front-illuminated chips. The differences are less than 40 eV for the 
back-illuminated chip, and although this is less than the ACIS back-illuminated spectral 
resolution of 100 - 150 eV, the CTI correction provides a significant improvement even for 
this chip. 

The issue of background was not straightforward for this observation since the 
southeast corner of the supernova remnant nearly fills the S3 chip. Using background from 
an off-source area and from blank-sky files did not change the fit results (see Table 1) 

1 Or visit the web site at http://heasarc.gsfc.nasa.gov 

2 Available online at http://www.astro.psu.edu/users/townsley/cti/ 
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within the errors. This is partially due to the small contribution (less than 0.5%) of the 
background to the source counts for RCW 86. Since standard blank sky background is not 
compatible with the CTI-corrected spectra, we have chosen to take background spectra from 
off-source areas on the same chip. For the front-illuminated 13 chip with more substantial 
CTI effects, varying CT1 corrections for different regions of the chip constrained us to 
choose for background a part of the same CTI region as used for the signal, rows 641-768. 
While these regions likely contain source as well as background photons, since they are from 
the same chip and have been processed identically with the source regions they should not 
introduce large errors into the source spectra. Residual instrumental background features 
such as a faint Si Ka line are present in at least one of our background-subtracted spectra, 
but since they are very weak they do not pose any problems for spectral fitting. 

3. RESULTS 
3.1. Images 

Figure 1 shows the image from the 12, 13, and S0-S3 chips of the SW corner of RCW 86, 
binned to a resolution of 2" and including the full energy range. The extensive filamentary 
structure is quite striking and persists down to the scale of the resolution. Most of this 
structure is parallel to the presumed edge of the remnant. However, some filaments can be 
seen that are almost radial in orientation, and some emission is considerably more diffuse. 
Compared to the ROSAT HRI image (Bocchino et al. 2000; Dickel et al. 2001), Figure 1 
shows emission extending considerably further ahead of the bright X-ray filaments in the S 
and SW. Some of that emission is bounded by a sharp-edged filament, but some (in the 
extreme SW corner, for instance), is not. 

Figure 2 contrasts the X-rays in soft (0.5 - 1 keV; red) medium (1-2 keV, green), and 
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hard (2-8 keV, blue) energy bands. (The soft and hard band images are shown separately 
in Figures 3 and 4, respectively.) The brightest emission is clearly relatively soft with 
some even softer emission further from the remnant center. The network of soft filaments 
is remarkably similar to the Ha image (Fig. 5, Smith 1997). The hard X-rays have a 
dramatically different morphology, with completely orthogonal filaments in some locations 
on S3 and more diffuse emission in 13. The brightest region in the hard X-ray band is the 
somewhat localized knot near the center of the 13 chip. These morphological differences 
strongly suggest a different origin for the hard X-rays than for the soft and medium bands, 
a conclusion confirmed by spectral analysis (see below). 

Radio and X-ray morphologies are compared in Figure 6, where X-rays below 2 keV 
are indicated with green, above 2 keV with blue, as in Figure 2, and radio (the Australia 
Telescope image: Dickel et al. 2001) in red. While neither X-ray band shows a structure 
closely resembling radio, there are clear correlations between the hard X-rays and radio: 
hard X-rays often appear at the edges of radio-emitting regions, most strikingly in a N-S 
filament of hard X-rays with radio extending beyond to the east, in the center of the S3 
chip. We have extracted radio and X-ray profiles from the three rectangular regions shown 
in Fig. 7, where we have averaged over the narrow dimensions of the rectangles. The profiles 
are compared in Fig. 8. All three regions show that thin X-ray filaments are associated 
with rises (or drops) in radio emission, including the hard knot in the center of 13. We 
interpret these as examples of shocks seen mostly edge-on, with the downstream direction 
the one in which radio emission rises. We shall argue below that these coincidences support 
a synchrotron origin for the hard X-ray emission. 
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3.2. Spectra 

Figure 7 shows extraction regions for the nonradiative Balmer filament, for the entire 
hard region on the S3 chip, for the synchrotron filaments alone, and for one region on the 13 
chip. The spectrum of the nonradiative shock, visible in optical as a bright Balmer filament, 
is soft (Figure 9), without any measurable emission above 4 keV. In particular, there is 
no Fe Ka line emission associated with this shock. We modeled its spectrum with plane 
shock models pshock and vpshock which are available in XSPEC (these models, described 
by Borkowski, Lyerly, & Reynolds 2001, are based on updated thermal plasma code of 
Hamilton & Sarazin 1984a). Because of calibration problems for the ACIS S3 detector at 
photon energies below 0.5 keV, we ignored data below this photon energy threshold in all 
our fits. 

We also excluded channels with energies from 1.17 keV to 1.31 keV in our fits to the 
X-ray spectrum of the soft region. This is the spectral region with "missing" Fe-L shell lines 
(Brickhouse et al. 2000). Modeling of the Fe L-shell complex in the equilibrium-ionization 
MEKAL and all non-equilibrium ionization (NEI) models in XSPEC is based on calculations 
by Liedahl et al. (1995), which did not include electron excitations to levels higher than 
n = 5. Brickhouse et al. (2000) found that Fe L-shell lines originating at these highly 
excited levels cannot be neglected in the ASCA spectrum of Capella. In Capella these lines 
are produced by 6 x 10 6 K (0.5 keV) plasma in collisional ionization equilibrium, and they 
cluster near 1.2 keV. Because plasma temperatures are similar in Capella and in our soft 
regions associated with Balmer filaments, we excluded channels near this energy where we 
expect Fe L-shell lines originating at n > 5 which are missing in the NEI spectral code. 

Results of plane shock fits to the soft spectrum are tabulated in Table 1. In the first fit 
with wabs + pshock models, relative abundances of heavy elements are fixed at solar values, 
but their absolute abundance with respect to H and He is allowed to vary. In the second 
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fit with wabs + vpshock models, abundances of individual elements are allowed to vary 
independently with the following exceptions: the C abundance is tied to the N abundance, 
and Ca and Ni are tied to Fe. The second fit reproduces the data reasonably well (reduced 
X 2 = 2.7, see Figure 9). Both models underpredict emission at energies below 0.5 keV by 
up to a factor of 2, which may be either caused by calibration/instrumental problems or 
by the presence of additional gas with much lower temperatures in the spectral extraction 
region. Electron temperatures are relatively low, 0.41 - 0.43 keV, shock ionization ages are 
relatively short, (3 — 5) x 10 10 cm -3 s, and elemental abundances are close to solar. The 
largest departures from solar abundances occur for Fe and Mg, whose apparent subsolar 
abundances might be caused by moderate depletion onto dust grains and by overlapping Fe 
L-shell lines from high n levels, which are adjacent to the Mg Ka line. An apparent excess 
of emission near 3 keV hints at the presence of the Ar Ka line (Ar is not included in the 
NEI code). 

There is no evidence for the presence of excess emission at high energies, suggesting 
that thermal plane shock models with constant electron temperature and approximately 
solar abundances provide a satisfactory description of the soft spectrum. This finding 
fully confirms our previous inferences (BRRD) about the thermal nature of the soft X-ray 
emitting filaments which are associated with the blast wave propagating into the ambient 
ISM. 

We realize that our very simple single plane-shock model for the soft emission is 
quite unrealistic. In the multiband X-ray image (Figure 2), it is clear that the bright 
thermal filaments have spectral differences on spatial scales much smaller than the size 
of our spectral extraction regions, and extending down nearly to the Chandra spatial 
resolution limit. The small values of ionization timescale r we find may be artifacts of 
using constant-temperature shock models to fit regions with clearly varying temperatures. 
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The focus of the present investigation is understanding the harder emission, and for that 
purpose we simply require an adequate description, if not a full explanation, of the soft 
emission. We defer to a later work a more extensive exploration of the nature of the soft 
emission and possible explanations for the short ionization timescale. 

Spectra of hard regions are shown in Figure 10 (the entire synchrotron-emitting 
region in the NE corner of the S3 chip) and Figure 11 (synchrotron filaments on S3 and 
13 chips). (See Figure 7 for the locations of these regions.) Results of fits with the srcut 
+ low-temperature vpshock models are shown in these Figures and tabulated in Table 1. 
The Fe Ka line is modeled by an additional pure Fe vpshock model (see § 4.2 for further 
details). Ideally, one uses the srcut model with radio information as input: the 1 GHz 
radio flux of the region and the spectral index. The integrated spectral index of the entire 
RCW 86 remnant is about 0.6, based on single-dish observations (typical errors are ±0.05), 
but accurate fluxes cannot be determined because the only available radio images are 
interferometer images, in which the absence of very short antenna baselines means that very 
smoothly distributed flux is resolved out. In fact, the MOST image at 843 MHz (Whiteoak 
& Green 1996) contains less than 40% of the single-dish flux at that frequency, while the 
ATCA image (Dickel et al. 2001) contains about 70% of the 1.3 GHz single-dish flux. In 
local regions, the total flux can be in error by even larger factors. As a result, we cannot 
use the measured radio fluxes to constrain the X-ray modeling. Instead, we have fixed the 
radio spectral index at 0.6 and actually fit for the radio flux. The values we obtain are quite 
reasonable and within factors of 2 of the fluxes measured from the radio maps (which differ 
between themselves by a comparable factor). 

Nonthermal X-ray emission from RCW 86 was reported from Rossi X-ray Timing 
Explorer (RXTE) observations by Allen et al. (1999). The spatial resolution of RXTE 
is only about one degree, so it cannot resolve RCW 86. Allen et al. (1999) report that 
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its spectrum observed with the PCA instrument is well described by a power-law with 
energy index T = 3.3 ± 0.2 between 10 and 24 keV. Our srcut model, with parameters fit 
to the hard region as described above, is slightly curved in that interval, with a mean slope 
between 10 and 20 keV of T = 3.5. It predicts a flux at 10 keV smaller than that quoted 
by Allen et al. (1999) by a factor of 12.5, and gives a 1 GHz radio flux of 2.59 Jy. If we 
simply scale up our X-ray model to the flux level reported by Allen et al. (1999) by that 
factor, the predicted radio flux at 1 GHz is 32 Jy. Green (2001 3 ) reports the observed 1 
GHz flux of RCW 86 to be 49 Jy. We regard this crude agreement, to within 50% after 
an extrapolation of nine orders of magnitude in photon frequency, as evidence of broad 
consistency between the whole-remnant spectrum seen by RXTE and the local spectra we 
fit in individual synchrotron filaments, consistency which represents strong support for the 
idea of a substantial synchrotron component in the spectrum of RCW 86. 

For the soft thermal regions, we also attempted to bound any nonthermal emission. 
We obtained a quantitative upper limit on the presence of a synchrotron component by 
adding to the fit an srcut component, with a 1 GHz flux determined from the Australia 
Telescope observations (Dickel et al. 2001) and assuming a spectral index of 0.6, allowing 
folios to vary. The fits rejected values of z/ ro iiofr larger than about 2 x 10 16 Hz. 

Attempts to fit the hard continuum and the Fe Ka line with pure thermal models 
are discussed below; we find that the resulting extremely short ionization timescales are 
inconsistent with the known age of the remnant, which should be the shock age if the Fe 
Ka line is due to solar abundance Fe from the ISM. 



3 Available online at http://www.mrao.cam.ac.uk/surveys/snrs/ 
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4. DISCUSSION 

4.1. Morphological Evidence for Synchrotron Emission 

We believe that the striking difference in morphology of hard and soft X-rays is 
strong evidence in favor of the synchrotron interpretation of the bulk of the harder X-ray 
continuum (above 2 keV). The spectrum of regions bright in X-rays below 2 keV is well 
described by thermal plane-shock models. The hard X-rays show a dramatically different 
structure, often occurring on the edges of regions of radio emission. It is sometimes claimed 
that a nonthermal interpretation of X-rays requires that the X-ray and radio morphologies 
be identical. This is only true to the extent that synchrotron losses on electrons can be 
neglected, rarely the case except for very young remnants such as the historical Galactic 
remnants. In a magnetic field of 10 /zGauss, electrons radiating the peak of their synchrotron 
emission at hu = 4 keV lose half their energy in 900 years, far less than the presumed age 
of order 10 4 yr, while electrons radiating chiefly at 1 GHz will last for over 10 7 yr. Thus 
we should expect synchrotron X-rays to be seen near the shocks where the electrons are 
presumably accelerated, while radio emission can be seen there as well as far downstream. 

In the complex inhomogeneous environment of RCW 86, shock structures can be 
found in various locations inside the outermost X-ray and radio emission. Some may be 
parts of the reverse shock or secondary shocks driven into denser clouds, and projection 
effects further complicate the interpretation. However, a particularly clear example of what 
appears to be an internal shock seen edge on is indicated in the central ("hard-filament") 
extraction region on chip S3 in Figure 7, and profiles along the indicated direction are shown 
in Figure 8. Here, both radio and X-ray emission appear fairly abruptly; the X-ray emitting 
region is only a few tens of arcseconds wide, while radio emission persists much further, 
just as expected if this is an outward facing shock at which electrons are newly accelerated. 
At a distance of 3 kpc, 1" corresponds to 0.015 pc or 4.5 x 10 16 cm. Shock velocities from 
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optical or X-ray analysis are of order 600-900 km s" 1 , so that one might expect post-shock 
flow velocities in the plane of the sky to be of the same order, but somewhat smaller. 
In 900 years, a flow at 500 km s _1 travels a distance / = 0.46 pc corresponding to an 
angular distance of 32". This compares reasonably well with the thickness of the overall 
X-ray emitting peaks in the profiles of Figure 8, allowing for projection effects and effects 
of magnetic-field geometry, which can account for the irregular shapes of the profiles. In 
general, we should expect synchrotron-emitting regions to have thicknesses of order a few 
tens of arcseconds, though if a shock is seen close to face-on, one will obtain more diffuse, 
fainter extended emission. These properties do appear to characterize the structure of the 
hard X-ray emission as well as can be expected. 

We do not see clear evidence for hard X-ray emission from the extreme outer edge of 
the X-ray emission seen in the lower part of the S3 chip, where at least one part of the blast 
wave is presumably encountering undisturbed upstream material. However, the overall 
faintness of the emission (presumably mostly thermal) is such that a faint nonthermal 
component might be hard to discern. Unfortunately, the emission there is too faint for 
useful spectral analysis with these data. Interior shocks could be portions of the outer blast 
wave encountering denser material and seen in projection, or portions of a reverse shock 
moving back into the interior from such an encounter. It seems that such interior shocks are 
where most hard X-rays originate, though not at the Ha-bright nonradiative shocks which 
are probably the blast wave moving into partially neutral material. This is an extremely 
interesting result, as it suggests that electrons may be accelerated at reverse shocks as well 
as at forward shocks - perhaps preferentially at reverse shocks in RCW 86. In the prototype 
synchrotron- X-ray-dominated remnant, SN 1006, no obvious reverse-shock structure is seen, 
and the hardest emission appears to come from closest to the outer remnant edge. While a 
reverse shock is almost certainly present in SN 1006, one can explain its absence in radio 
or X-ray nonthermal emission if the magnetic field threading the ejecta is very weak, as 
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one might expect if it is just the magnetic field of the progenitor, enormously diluted by 
expansion with flux freezing. The situation appears to be different in RCW 86. 

4.2. Spectroscopic Evidence for Fe-rich SN ejecta 

The Fe Ka line at 6.4 keV in RCW 86 has been puzzling since Vink et al. (1997) 
reliably determined its energy, because such a low energy implies the presence of only very 
low ionization stages of Fe in the X-ray emitting plasma. Vink et al. (1997) suggested strong 
deviations from a Maxwellian electron distribution in order to account for both the presence 
of the Fe line and the weakness of other emission lines at low energies. Subsequent studies 
(BRRD; Bamba, Koyama, & Tomida 2000; Bocchino et al. 2000) have generally postulated 
its origin in a hot (~ 5 keV) plasma in the remnant's interior. A substantial Fe Ka strength 
suggests that Fe is overabundant by a factor of several with respect to cosmic (solar) 
abundances, but this finding is somewhat uncertain in view of the poor photon statistics 
and uncertainties in atomic data. Although our Chandra observations of RCW86 also suffer 
from poor photon statistics at high energies, the high spatial resolution of Chandra allows 
us to spatially separate the synchrotron-emitting regions and the soft thermally-emitting 
X-ray filaments. We find that Fe Ka is associated with the synchrotron-emitting regions, 
and not with the thermal filaments. Apparently, fast electrons capable of exciting Fe Ka 
are present in synchrotron emitting regions, and absent in soft X-ray filaments. 

It seems natural to assume that the Fe Ka line is produced in a hot thermal plasma 
with an approximately cosmic composition, where most free electrons are provided by the 
ionization of H and He. In this hypothesis the observed high-energy continuum must be 
thermal, because a substantial contribution from synchrotron emission would demand a 
high derived Fe abundance for the line to rise above the continuum, and this would be 
obviously in conflict with the assumed solar composition. This interpretation unexpectedly 
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encounters a timescale problem which makes it untenable, as we now demonstrate. We use 
the fitted ionization timescale n e t, and derive n e from the observed X-ray flux in the model 
to obtain a characteristic age of the shocked plasma under this assumption. 

For a plasma consisting mostly of H and He, we can set a lower limit to the electron 
density in the synchrotron-emitting region by fitting a thermal NEI model to the high-energy 
Chandra spectrum. But because the Chandra spectrum is rather noisy in the vicinity of the 
Fe Ka line, we shall obtain the temperature, Fe abundance, and ionization timescale from 
high-quality ASCA GIS spectra presented by us previously (BRRD). (The analysis of the 
Fe Ka line in BRRD was influenced by an error in all NEI models in XSPEC vlO, which 
resulted in a missing Fe Ka fluorescence line at extremely low ionization timescales relevant 
for RCW86. This error is not present in the NEI models in XSPEC vll.) The corrected 
values for the vnei model fit to GIS spectra above 5.5 keV are: plasma temperature of 5.0 
keV (with a 90% confidence interval of 3.5-7.6 keV), a very low ionization timescale r of 
2.8 x 10 8 cm -3 s (less than 6.6 x 10 9 cm -3 s with 90% confidence), and an Fe abundance 
of 1.6 solar (with 90% confidence interval 1.1-2.5). By using these best-fit values, we 
find that this model must be normalized to 0.0048 (in XSPEC units for thermal models: 
normalization = (10" 14 / 'An D 2 ) J n 2 dV, where D is the distance to the source in cm and n is 
the electron density). This value is required to account for the high energy emission seen in 
Figure 10, implying an emission measure EM = J n e nudV of 5.2 x lO 56 ^ cm" 3 (where d% is 
the remnant's distance in units of 3 kpc). The extraction region is approximately 2(4 x 5f0 
in size, spanning the entire width of the synchrotron-emitting region. This translates into 
a surface area of 2.1 x 4.2 pc 2 at 3 kpc distance. Assuming 4.2 pc depth, we estimate the 
volume of the X-ray emitting region to be 37c?| pc 3 . For H and He dominated plasma, we 
then arrive at an electron density n e of 0.76 f^^d^ cm" 3 , where / is the volume filling 
fraction of the hot gas. The plasma characteristic "age" is then < r/n e = 12 f l ^ 2 d l J 2 yr, 
which is far too short compared with the age of the remnant. 
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The upper limit on r from fits to ASCA GIS data is of course much larger than 
the best-fit value of 2.8 x 10 8 cm -3 s used in these estimates, but we can establish a 
more stringent upper limit on r which applies for plasmas with approximately cosmic 
composition. We noted previously an extreme sensitivity of NEI models to ionization 
timescale (Figure 8 in BRRD). At extremely low ionization timescales, most abundant 
elements such as O have not been stripped of their L-shell electrons, resulting in weak 
X-ray lines in the Chandra ACIS spectral range and low X-ray fluxes at low energies. As 
the ionization timescale increases, electrons are stripped from their L shells, and the soft 
X-ray flux increases dramatically. By varying the ionization timescale, we find that the vnei 
model with the temperature and Fe abundance quoted above and Nh = 4.2 x 10 21 cm~ 2 
(from model B in Table 1) produces strong O lines at low energies which are inconsistent 
with Chandra spectra for r > 1 x 10 9 cm' 3 s. (We obtained this conservative upper limit by 
requiring that the calculated spectrum not exceed the observed count rate at the location of 
O lines. But in contrast to observations, O lines are very prominent in strongly-underionized 
models, so based on consideration of the spectral shape one can probably reject models 
with r as short as 5 x 10 8 cm~ 3 s.) This constraint on ionization timescale gives an age 
of less than 40 yr. By assuming a much higher plasma temperature we can perhaps push 
this limit to ~ 100 yr, which is still much less than expected for a relatively old remnant 
such as RCW 86. This timescale problem is not an entirely new finding, as an inspection of 
previous results obtained with various satellites demonstrates it clearly, but new constraints 
from Chandra are particularly severe. The Fe Ka line and continuum at high energies 
cannot be easily produced by a thermal plasma with approximately solar abundances. 

A consideration of Coulomb interactions between ions and electrons (see the Appendix) 
leads us to a conclusion that the observed continuum in regions of the remnant with hard 
spectra cannot be produced by nonthermal bremsstrahlung in cosmic abundance plasmas. 
This leaves synchrotron radiation as the preferred explanation for the observed continuum 
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in the whole Chandra spectral range, including high energies. But the Fe Ka line must be 
produced by collisional excitations and ionizations by much less energetic thermal electrons. 
Because the Fe Ka line strength relative to the high energy continuum already implies 
above solar abundances, it is likely that we are observing highly-enriched SN ejecta. SN 
ejecta were already considered in some detail by Bocchino et al. (2000), but now we have 
indirect evidence for the presence of very Fe-rich ejecta. For simplicity, we considered 
pure Fe-rich ejecta and modeled them by the vpshock model with temperature 5 keV and 
ionization timescale 10 9 cm -3 s. These plasma parameters are just rough guesses: at much 
lower temperatures most thermal electrons are not capable of exciting the Fe Ka line, while 
at ionization ages above a few xlO 9 cm -3 s, Fe L shell emission becomes very prominent 
(in conflict with Chandra observations). In addition to the pure Fe vpshock model which 
accounts for the Fe Ka line, we used the srcut model for modeling synchrotron emission. 
The blast wave seen in projection against synchrotron-dominated regions is modeled 
by a low-temperature vpshock model, either with solar abundances, or with abundances 
determined from our fit to the spectrum of the Balmer-dominated nonradiative shock. 
Results are tabulated in Table 1. 

The Fe Ka line is unambiguously detected only in the spectrum of the entire 
synchrotron-emitting region on the NE corner of the S3 chip (see Figure 7). The emission 
measure EM Fe = J n e np e dV of pure Fe-rich ejecta is 3.1 x 10 52 <i 2 cm" 3 , based on the 
normalization of model B in Table 1. With the estimated X-ray emitting volume V of 37rfg 
pc 3 , we arrive at n e = 0.017(n e /10n Fe ) 1 / 2 /~ 1/2 ^3 1/2 cm" 3 . The factor (n e /10n Fe ) 1/2 is of 
the order of unity for a pure Fe plasma with r = 10 9 cm" 3 s, but it is somewhat uncertain 
because of the poorly known electron distribution function and a possible contribution to n e 
from other heavy elements which might be present in the ejecta. The electron distribution 
function may well be non-Maxwellian in such a low-density and high-temperature plasma 
because of the long Coulomb collision timescale. In particular, the distribution may be 
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bi-modal, with free preshock electrons heated to high temperatures at the collisionless 
shock, and with much colder electrons released in the process of collisional ionization in the 
postshock region (Hamilton & Sarazin 1984b). The electron density n e just quoted refers 
to the hot component, and not to cold electrons with energies much below the Fe K-shell 
ionization threshold. For a pure Fe plasma with electron density 0.017 cm -3 and r = 10 9 
cm -3 s, the plasma age is about 2000 yr, which is consistent with the relatively large age of 
RCW 86. 

The presence of low-density Fe-rich ejecta in RCW 86 is consistent with the large 
observed equivalent widths of the Fe Ka line, which already led Bocchino et al. (2000) 
to suggest the possible presence of SN ejecta in this remnant. The morphology of 
synchrotron-emitting filaments apparently associated with the ejecta is also suggestive. The 
X-ray and radio synchrotron-emitting filaments are located mostly interior to the bright 
nonradiative and radiative shocks in the "knee" region of the remnant. These filaments form 
a ridge which extends much further to the east, and which is clearly located well within the 
blast wave outlined by nonradiative Balmer-emitting shocks (Smith 1997). The brightest 
part of this ridge is in the region of highest pressures where the most intense interaction of 
supernova ejecta with the ambient ISM is occurring, as judged by the presence of bright 
radiative and nonradiative shocks. This morphology is consistent with identification of this 
ridge with the reverse shock. In addition to solving the short timescale puzzle, the presence 
of low-density Fe-rich SN ejecta also removes a need for a mixed thermal-nonthermal 
continuum postulated by us previously for the synchrotron emission region (BRRD; see 
also Bamba, Koyama, & Tomida 2000). This required a coincidental match in intensities 
of synchrotron and thermal continuum emission in the 1-10 keV range. This is no longer 
necessary, since the observed continuum is mostly synchrotron emission, and the thermal 
contribution to the continuum may be as much as 2 or 3 orders of magnitude fainter. 
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4.3. X-Ray Synchrotron Spectra 



No synchrotron component was required to obtain a reasonable fit for the soft region 
(Fig. 9); an upper limit to the rolloff frequency of about 2 x 10 16 Hz was obtained as 
described in the previous section. We find values of f ro iioff of (7 — 10) x 10 16 Hz from the 
fits to hard regions, with remarkably little dispersion. These values are two to three times 
higher than those we reported in BRRD, for two reasons. First, the enormously improved 
spatial resolution of Chandra allows us to select spectral regions for extraction that are 
localized to the very hardest X-ray emission, whereas the fits based on ASCA data are 
necessarily averages over much larger regions, lowering the mean value of v ro \\ $. Second, as 
described above, we now believe the entire hard continuum, rather than just part of it, to 
be due to synchrotron X-rays. 

The values of z/ ro iioff can be understood in the framework of simple estimates of particle 
shock-acceleration times (see, for example, Reynolds 1998). If the maximum energy that 
electrons can reach is limited by radiative losses, the energy at which energy gains by 
acceleration equal losses due to radiation is about 



where 77 is the ratio of electron mean free path to gyroradius (and is assumed constant), B 1 
is the upstream magnetic-field strength, and u 8 is the shock speed in units of 10 s cm s _1 . 
We have assumed a shock compression r of 4. The factor Rj is Jokipii's (1987) expression 
for the obliquity-dependence of the particle acceleration rate; where the shock is nearly 
perpendicular (i. e., the shock normal nearly perpendicular to the upstream magnetic field 
direction), 



so that for 77 substantially greater than 1, 77 Rj oc rj 1 and E m oc r] +1 / 2 . Electrons of this 



E m ~ 50 (vRj)- l/2 (5i/3 /iG)" i/2 M8 ergs, 
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energy radiating in the downstream field B 2 produce their peak power at a frequency 

z/ c ~5xl0 16 (r ] R J y 1 (B 2 /4 B l )u s 2 . 

If reverse shock velocities are comparable to the blast-wave velocities of 400 - 900 km s -1 , 
then, we can reach the observed values of z/ rolloff for compression ratios substantially larger 
than 4, and/or values of rj substantially greater than 1. (Notice that the characteristic 
frequency radiated by electrons whose maximum energy is set by radiative losses is 
independent of the magnetic-field strength.) It should not be surprising that the hard 
emission is so restricted in spatial distribution, then; many of the shocks present in RC W 
86 will not be fast enough to produce synchrotron X-ray-emitting electrons. Only in a few 
regions where the conditions are right will we expect such electrons to be accelerated. 

We should add that the model we use to describe the synchrotron emission, srcut, 
is highly oversimplified: it is just the synchrotron emissivity of a power- law distribution 
of electrons with an exponential cutoff (the fastest plausible cutoff for shock-accelerated 
particles; see references in Reynolds & Keohane 1999), radiating in a uniform magnetic 
field. Much more elaborate models of spherical remnants in which radiative losses limit 
electron acceleration were described in Reynolds (1998); those models, not available in 
XSPEC, produced synchrotron spectra with a somewhat different shape than srcut. While 
the curvature within the Chandra bandpass is not enough to distinguish the difference, the 
inferred values of f ro iiofr would probably be slightly different for a true loss-limited model. In 
fact, since any spatial inhomogeneities will broaden the cutoff, srcut represents the steepest 
possible cutoff of emission, and will give the highest rolloff frequency compatible with the 
observed data. Any more realistic model should produce a somewhat lower value. 

Our fits to the synchrotron component could be substantially better constrained if 
radio fluxes were known more accurately. The most straightforward way to accomplish 
this would be to map RCW 86 with a large single dish, such as the Parkes telescope, and 
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combine those data with the ATCA interferometer image. This would allow considerably- 
firmer determinations of f ro iiofr- 

We have argued that the hard continuum and Fe Ka line tend to be found together, 
strongly suggesting that it is the reverse shock that is responsible for accelerating particles 
to the highest energies. This conclusion is important to confirm, since it seems not to be 
obviously the case in other remnants with synchrotron X-ray emission. There is no obvious 
reason that reverse shocks should not be as effective as the blast wave at accelerating 
particles, for the same values of physical parameters such as shock speed and 77. In RCW 
86, blast-wave emission, while detectable, is so faint that synchrotron X-rays would be 
below the level of background, so we cannot rule out electron acceleration to comparably 
high energies at the blast wave. 

Observing particle acceleration at the reverse shock may also have interesting 
consequences for the acceleration of cosmic rays. The direct acceleration of enriched 
material in supernova remnants might play a role in models of Galactic cosmic-ray 
acceleration and of chemical enrichment of the light elements Li, Be, B in the early Galaxy, 
when the carbon and oxygen from which these elements are spalled were less abundant in 
the general interstellar medium. 

Even more intriguing is the possibility that the reason RCW 86 shows particular 
evidence for particle acceleration at the reverse shock is that the radiating particles are not 
electrons but positrons produced by decay of radioactive 56 Ni synthesized in the explosion. 
A longstanding theoretical difficulty in electron acceleration has been the "injection 
problem" : thermal electrons have such small gyroradii that they see the shock as a smooth 
transition and not the discontinuity in flow speeds required for standard diffusive shock 
acceleration (see, e. g., Levinson 1992, 1994 for extensive discussion). Various suggestions 
have been made about how Nature evidently overcomes this difficulty. However, it has been 
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noted (Ramaty & Lingenfelter 1979; Ellison, Jones, & Ramaty 1989) that positrons released 
by radioactive decay of 56 Ni and 56 Co have typical energies of order 1 MeV, where they 
should see the shock as a discontinuity and be accelerated as efficiently as protons of similar 
energies. A reverse shock moving into Fe-rich ejecta should encounter a population of 
"pre-injected" energetic positrons which can be immediately accelerated efficiently. While 
this is a possibility difficult to rule out, we note that the decay lifetimes of 56 Ni and 56 Co 
are both less than a year; such positrons were almost entirely produced about 10 4 years 
ago, just after the supernova. Any X-ray-synchrotron-emitting positrons cannot have been 
accelerated soon after that, because the radiative lifetime is less than 1000 yr. However, 1 
MeV positrons will lose energy primarily by adiabatic expansion losses, which could reduce 
them to no better than thermal electrons after 10 4 yr. This fascinating possibility must 
remain a speculation at this time. 

5. CONCLUSIONS 

The major conclusion of our paper is that the detailed morphology of soft and hard 
X-rays in RCW 86 strikingly supports the hypothesis of different origins. Our spectral 
analysis confirms that the hard X-rays are best described as a synchrotron continuum 
with Fe Ka emission that must come from a thermal component with a substantial 
overabundance of iron with respect to solar. This result implies that some of the original 
supernova ejecta are still unmixed after ~ 10 4 yr, and that the reverse shock into those 
ejecta can also accelerate electrons to X-ray synchrotron-emitting energies of order 50 TeV. 

Both these conclusions substantially complicate the task of understanding the X-ray 
spectra of middle-aged supernova remnants. Apparently even after ages of order 10 4 years, 
remnants are not in the simple Sedov phase, and the possibility of synchrotron X-ray 
emission confusing the analysis of thermal emission cannot be neglected. The only ultimate 
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solution to this problem requires X-ray spectroscopy with much higher energy resolution. 
(Grating resolution is adequate, but grating observations of large extended sources with 
Chandra or XMM-Newton are virtually uninterpretable.) The Astro-E2 microcalorimeters 
hold out the possibility of spatially resolved spectroscopy with enough resolution to allow 
the inference of temperatures from line complexes alone without requiring any assumptions 
about the continuum. In that case, one can immediately determine the extent of synchrotron 
contributions to the continuum. 

With Chandra's spatial resolution, we can for the first time discern discrete shock 
features and minimize projection effects, which should aid substantially in the detailed 
thermal modeling that will be the next step in understanding RCW 86. The hard X-ray 
emission appears to delineate particular shocks, where both radio and X-ray emitting 
electrons are accelerated. The greater extent of radio emission behind such features 
is expected from preferential energy losses on the much more energetic X-ray emitting 
electrons. The width of the X-ray emission is quite consistent with this interpretation, for 
expected post-shock velocities of hundreds of km s -1 and magnetic field strengths of order 
10 /i-Gauss. 

The observed rolloff frequencies in areas containing synchrotron continua can be 
understood in terms of radiative-loss-limited shock acceleration, but simple estimates 
require favorable conditions to reach the inferred electron energies. In particular, nearly 
perpendicular shocks reach higher energies more readily, and this may be one factor that 
distinguishes shocks that produce synchrotron X-ray continua from those which do not. 
We repeat, however, that the srcut model is highly simplified, and fits with more elaborate 
models might produce somewhat different values for ^ ro ii ff- 

Progress in understanding RCW 86, and the general problems it raises for the 
understanding of supernova remnants, requires further observations in various bands. 
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XMM-Newton observations can obtain significantly higher signal-to-noise ratios at high 
energies, allowing testing in other regions of our claim of an association of Fe Ka emission 
with regions of hard continuum. Optical and infrared observations in principle should 
be able to detect synchrotron continuum at a level exactly predictable from our fits, but 
distinguishing it from the confusing emission might prove extremely challenging. Since our 
fitted values of z^oiioff are far above optical frequencies, we expect that there should be little 
difference between the optical or IR morphology of synchrotron emission and the radio 
morphology. A strategy to search for extremely faint diffuse emission then might center on 
sharp features in the radio image, such as the edge coincident with a hard X-ray filament in 
the middle of chip S3. 
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A. Appendix 

Here we demonstrate that nonthermal bremsstrahlung is unlikely to be important in 
the regions we have analyzed in RCW 86. First, we shall show that such an explanation for 
the lower-energy continuum would require a huge unseen cold-electron population of which 
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the observed ~ 2 keV electrons are the tail. The short timescale for Coulomb equilibration 
of these ~ 2 keV electrons with the cold population would require that the cold population 
have an enormous density; it would itself cool rapidly and produce copious unseen optical 
emission, and the energy required to maintain the warmer population would be exorbitant. 
We then show that even the ~ 5 keV continuum near the Fe Ka line is unlikely to be 
nonthermal bremsstrahlung. Again, the energetics of maintaining the hot population 
against energy transfer into the cold population would require that a significant portion of 
the blast energy go into electrons with energies of order 5 keV, which seems unreasonable. 

We showed in §4.2 that the Fe Ka line cannot be produced in a thermal plasma with 
approximately solar abundances. Is it possible to account for the observed Fe Ka emission 
by postulating deviations from a Maxwellian electron distribution in such plasma, as 
suggested by Vink et al. (1997, 2002)? First we consider whether the featureless spectrum at 
lower energies (~ 2 keV), which we interpret as nonthermal synchrotron emission, could be 
attributed to nonthermal bremsstrahlung produced by suprathermal electrons accelerated 
in the shock wave. The following arguments based on order-of-magnitude estimates 
demonstrate that this is not possible. For the nonthermal bremsstrahlung hypothesis, one 
needs a large thermal pool of cold electrons with temperature T c , and much less numerous 
suprathermal electrons forming a power-law tail at high electron energies. The description 
of these fast electrons in terms of a power law is a crude approximation because of Coulomb 
losses to the cold electrons, and because the particle acceleration process at the shock may 
produce a curved (concave) spectrum instead of a power law (Ellison & Reynolds 1991; 
Reynolds & Ellison 1992). In view of these complexities, we do not attempt to model 
the electron distribution function, but assume instead a Maxwellian distribution of hot 
electrons. Such a two-component (cold+hot) thermal model is sufficient for our present 
purpose of estimating the number density of fast electrons which would be required to 
account for the observed strength of the featureless X-ray continuum. The bulk plasma 
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temperature T c must be low enough that the ionization state of the gas is so low that 
practically no X-ray lines are produced in the Chandra energy range, to be consistent with 
the observed featureless spectra. The observed X-ray continuum would then be free-free 
emission produced by fast electrons in collisions with protons and a particles. In order to 
estimate the required density of fast electrons, we fit a standard free-free emission model to 
the spectrum shown in Figure 10. We obtain a hot electron temperature of 1.45 keV, 
and corresponding emission measure EM h = J n e hriHdV of 2.6 x 10 57 g?§ cm" 3 (where n e h is 
the number density of hot electrons, which must be much smaller than the hydrogen density 
n H ). For the X-ray emitting volume of 37d 3 pc 3 estimated in §4.2, n e h = 2A/ f h n H dz cm -3 , 
where fh is the unknown volume filling fraction of fast electrons. 

The hot electron temperature of 1.45 keV implies an average energy E for fast electrons 
of 2.2 keV. Electrons with such a low energy lose it rapidly in Coulomb collisions with 
cold electrons in the thermal pool. The Coulomb energy loss timescale t E is equal to 
32£ 3 /7G((£/T c ) 1 / 2 )n e lnA yr (Spitzer 1962), where n e = 1.2n H is the density of cold 
electrons, In A is the Coulomb logarithm, and the function G is tabulated by Spitzer (1962). 
Because expected electron densities are high (n e ^ 1 cm -3 ), this is a short time scale. 
For example, for E = 2.2 keV and T c = 0.1 keV, we find In A = 30, (E/T c )^ 2 = 4.7, 
G(4.7) = 0.024, and t E = MOn" 1 yr. This short Coulomb energy exchange time t E implies 
that 2 keV electrons are strongly coupled to the low-temperature thermal plasma through 
Coulomb collisions. This strong thermal coupling demands that the thermal energy in fast 
electrons be much smaller than in the cold thermal pool, otherwise the energy so quickly 
transferred to the cold electrons would raise their temperature enough to produce strong 
X-ray lines, contrary to observations. 

Let R = (n e + nH)T c /n e hTh denote the cold/hot energy ratio, which we demand to 
be much larger than unity. Using our estimates for n e h obtained above from the strength 
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of the X-ray continuum, R = Q.§2f h n 2 e d?,T c /T h . Because of the strong coupling through 
Coulomb collisions, the volume filling fraction f h is expected to be low for shock-accelerated 
fast electrons, on the order of t^/t s or less, where t s is the shock age. The expected shock 
age t s should be a sizeable fraction of the remnant's age, which is ~ lO 4 ^ yr (Rosado et al. 
1996). Denoting t S)4 = t s /10 A yr and using the expression for obtained above, we arrive 
at R = 0.002E 1 / 2 T c n e /t sA G((E/T c ) 1 / 2 )\nA. By equating R to 1, we obtain a conservative 
lower limit to the cold electron density n e of 500t S)4 G(( J B/T c ) 1 / 2 ) \nkE~ l / 2 T~ l cm" 3 . This 
is a very large electron density, e. g., even for a relatively recent (t Sj4 = 0.2) shock with 
E = 2.2 keV, and T c in the range of 0.05-0.1 keV, we arrive at n e ~ 500 cm -3 . Even if 
T c < 0.1 keV, the density of this cold plasma is large enough to produce a significant amount 
of very soft X-ray emission which would be easily seen in X-ray spectra. In addition, such 
plasma is located close to the peak of the cooling function, and it would itself cool and 
recombine on timescales much smaller than the age of the remnant even in the presence of 
magnetic pressure support. Such a cooling and recombining plasma should produce strong 
optical emission. There is no evidence for optical emission at the location of the hard X-ray 
emitting filaments in the optical images presented by Smith (1997). Another problem is the 
high (6.5 x 10~ 8 dynes cm~ 2 ) pressure of hot electrons, which exceeds the estimated ram 
pressure of the blast wave by an order of magnitude or more. This problem becomes more 
acute if the filling fraction fh is not small. In this case one must assume that an additional 
energy source exists throughout the X-ray emitting plasma, which could heat hot electrons 
and balance their Coulomb losses. The power required to achieve this balance far outstrips 
the amount of kinetic energy in the blast wave. We conclude that the observed X-ray 
continuum emission at low energies cannot be produced by nonthermal bremsstrahlung. 

While the low-energy continuum is almost certainly synchrotron radiation, the 
high-energy continuum in the vicinity of the Fe Ka line could in principle be produced 
by nonthermal bremsstrahlung by high energy suprathermal electrons. Because of the 
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high (~ 5 — 10 keV) electron energies required to produce this continuum, such electrons 
will be less strongly coupled to the cold electrons through Coulomb collisions than the 
~ 2 keV electrons considered above. The high energy continuum is nevertheless quite 
strong, implying a high density n e h for these fast electrons and nonnegligible thermal 
coupling between cold and fast electrons. We can repeat our previous arguments where 
now the subscript h refers to the ~ 5 — 10 keV electrons. Using our thermal fit to the 
high energy continuum and our estimates of the X-ray emitting volume (see §4.2), we 
arrive at n e h = 0.48/ fan nds cm -3 . Next, we estimate the energy transfer rate from hot to 
cold electrons, again approximating their electron distribution functions by Maxwellians 
with low (T c ) and high (7^) temperatures. This rate (per unit volume) is equal to 
— \n e hkdT h jdt = ^n e hk(T h — T c )/t eg , where t eq is the equipartition time between two 
Maxwellian electron distributions given by Spitzer (1962). With our estimate for n eh , 
and assuming that T c <C T^, we arrive at the energy transfer rate of 0.86kTh/n e t eq fhd3. 
Note that because t eq oc n" 1 , this energy transfer rate is independent of the density of 
cold electrons. For kT h = 5 keV, fa — 1, we arrive at a conservative lower limit to the 
energy transfer rate of 3.4 x lO -18 ^ 1 ergs cm' 3 s -1 . If we now neglect fast electron energy 
losses, i. e., if we assume that T h is constant so that the fast electrons coexist with the 
cold electrons throughout the whole X-ray emitting plasma, we find that in 1,000 yr the 
total amount of transferred energy is equal to 1.1 x 10~ 7 ergs cm' 3 . This would result in a 
thermal pressure of the cold plasma of 7.1 x 10~ 8 dynes cm -3 , more than the ram pressure 
of the blast wave, where velocities range up to 900 km s _1 and preshock densities may 
be as high as ~ 1 cm' 3 (Ghavamian et al. 2001). While there are some uncertainties in 
estimates of the preshock density because of the unknown preshock ionization fraction, this 
finding suggests that ~ 5 — 10 keV electrons would play a major role in the energy balance 
of the plasma. This is not expected in standard particle acceleration models where most 
of the shock kinetic energy is transfered to cosmic-ray particles, and would rule out the 
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simultaneous production of ~ 100 times as much energy in accelerated ions, as observed in 
cosmic rays at Earth. 

This problem remains if we relax our assumption of constant Th by allowing for 
Coulomb energy losses. We assume that fast electrons are injected at the shock front with 
energy E and then lose their energy on a timescale t E = 5.3 x lO 3 ^" 1 yr {t E was estimated 
here in the limit of T c = 0). For expected electron densities n e far in excess of 1 cm -3 , 
Coulomb energy losses are important. The volume filling fraction f\ for plasma containing 
fast electrons is then small, less than ~ £e/£ s = O.dSt'ln' 1 . With this volume fraction, we 
arrive at n eh = 1.6t s 4 cm -3 , and with T h = 5 keV the pressure of fast electrons is equal 
to 2.6 x 10~ 9 dynes cm" 2 even for a relatively young (t Sj4 = 0.2) shock. This pressure is 
substantial when compared with the estimated ram pressure, again implying a major role in 
the plasma energetics for ~ 5 — 10 keV electrons. It is hard to explain this in the framework 
of our current understanding of collisionless plasmas and the particle acceleration process. 
While we cannot completely rule out nonthermal bremsstrahlung as the explanation for 
the high-energy continuum on an observational basis alone, we do not favor it because of 
excessive energy requirements. 
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Table 1. Spectral Models 
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Fig. 1. — Total energy band (0.5-8 keV) Chandra image of RCW 86. The scale ranges from 
6 x 1CT 8 to 1 x 1CT 6 photons cm" 2 s" 1 arcsec -2 . The two chips in the left column are, from 
top (north) down, 12 and 13; in the right column, SO, SI, S2, and S3. This and following 
images have been binned to a resolution of 2". 

Fig. 2. — Mosaiced three-color Chandra images of RCW 86: red represents 0.5-1 keV pho- 
tons; green, 1-2 keV; and blue 2-8 keV. 

Fig. 3.— Soft (0.5-1 keV) X-ray image of RCW 86. 

Fig. 4.— Hard (2-10 keV) X-ray image of RCW 86. 

Fig. 5. — Ha image from Smith (1997) (in green) compared with the soft (0.5-1 keV) X-ray 
image (in red). 

Fig. 6.— The 1.4 GHz ATCA image (Dickel et al. 2001) (in red) compared with hard (2-10 
keV) and soft (0.5-1 keV) Chandra images (in blue and green, respectively). 
The radio image has been smoothed to a resolution of 8". 

Fig. 7. — Chandra image in the 0.5-8 keV energy band with spectral extraction regions and 
spatial profiles indicated (in red and yellow, respectively). 

Fig. 8. — Spatial profiles of hard (2-8 keV) X-ray emission and 1.4 GHz radio emission (solid 
and dashed lines, respectively) for three positions shown in Figure 7, from west (top) to east 
(bottom). X-ray emission is scaled independently for each profile relative to radio emission. 
Distance is measured in arcseconds, from E to W (top profile), and from NE to SW (middle 
and bottom profiles). 

Fig. 9. — ACIS S3 X-ray spectrum of the bright nonradiative Balmer filament. Plane-parallel 
shock model (model B from Table 1) is shown by solid line. 
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Fig. 10. — ACIS S3 X-ray spectrum of the entire synchrotron-emitting region. Model B 
from Table 1 is shown by solid line. Individual models components are also shown: srcut 
- short-dashed line, pure Fe vpshock model - long-dashed, low-temperature vpshock model - 
dash- dotted. 



Fig. 11. — ACIS S3 (top) and 13 (bottom) X-ray spectra of synchrotron-emitting filaments, 
with best-fit models from Table 1 shown by solid lines. Individual model components are 
also shown. 




Fig. 7 
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Fig. 8 
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